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o : 

^ We discuss influences of bulk viscosity on the Early Universe, which is modeled by 

Friedmann-Robertson- Walker metric and Einstein field equations. We assume that the 
matter filling the isotropic and homogeneous background is relativistic viscous character- 

Qized by ultra-relativistic equations of state deduced from recent lattice QCD simulations. 
We obtain a set of complicated differential equations, for which we suggest approximate 
solutions for Hubble parameter H. We find that finite viscosity in Eckart and Israel- 
Stewart fiuids would significantly modify our picture about the Early Universe. 
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1. Introduction 



$H i Bulk viscosity (BV) is of essential importance in Early Universe, because it can oc- 

cur in an isotropic and homogeneous matter. Although, BV likely vanishes at very 
high temperatures, it remains important in mixed phases of matter. It diverges near 
the phase transition(s). First-order theory for relativistic fluids has been suggested 
by Eckart^ and Landau and Lifshitz.^ A relativistic "second-order" theory has been 
introduced by Israel and Stewart (IS).^'^ Due to the complicated characters of re- 
sulting differential equations, very few solutions of gravitational field equations are 
known in framework of causal fluid theory. For technical reasons, most investigations 
for dissipative causal cosmologies are assuming symmetric Friedmann-Robertson- 
Walker (FRW) metric. 

. Recent RHIC results give a strong indication that hot, dense and viscous mat- 

ter has been produced^ in heavy-ion collisions. This matter likely agrees with the 
Lattice QCD simulations,^ which show that BV is not negligible near QCD critical 
temperature Tc.^'^ It has been shown that BV at high temperature T and weak 



' coupling as, £, ~ a^T^/lnaJ^.^ Following the prescription in Ref.,^ causal viscous 

fluid is given by ^ oc p, where p is energy density. 

In this letter, we assume finite bulk viscosity coefficient ^ in Eckert- and IS- 
fluids and apply QCD EoS to solve Freidmann's equation. We suggest numerical 
and approximate solutions for Hubble parameter H in dependence on comoving 
time t. We compare our results with the non-viscous fluid. Also We imply our 
assumptions to relate t with the temperature T. 

2. Model and Results 

Using natural units, k — c — 1 and /i = 27r, we assume that background geometry 
of the Early Universe is filled with bulk viscous fluid, which is given by a spatially 
flat FRW metric. In spherical coordinates, the line element 

ds^ = df - a^{t) [dr^ + {dd^ + sin^ 9d(l)^)] . (1) 
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At vanishing cosmological constant A, Einstein gravitational field equations in a flat 
Universe (fc — 0) read 



Rik — -^gikR — StT G Tik- 



(2) 



Inclusion of bulk viscous effects can be generalized through an effective pressure 
n, which is formally included in the effective thermodynamic pressure Peff,^ where 
Pe/ / = P + n and p is the thermodynamic pressure. Then, in comoving frame, the 
energy momentum tensor has the components Tq — p, = = Tg = —p^ff. For 
FRW line element, Eq. ([1]), field equations become 

3i?2 = p, 2H + 3H^ =peff, (3) 

where one dot denotes differentiation with respect to t. In a closed system, total 
energy density of cosmic matter is conserved, i.e., — and 

p + 3i/(pe//+p)-0, (4) 

where H — a/ a is the Hubble parameter. In presence of bulk viscous stress 11, 
Eq. (Ill) can be written as 

p + 3H{p + p) = -3nH. (5) 
According to Eckert theory, bulk viscous stress is given as 



n = -3S,H, 

which has in IS-theory,^''* following evolution: 
rri + n = -3£,H - irn I 3H 



(6) 



(7) 



where T is temperature and r is relaxation time. 

EoS for p and T can help to have a closed system from Eq. ([3]), (O and d?]). r 
and ^ can be determined according to phenomenological approaches. For instance, 
BV in quark-gluon plasma (QGP) at high T can be estimated from recent lattice 
QCD simulations.6'i° 

(3p^, ^ = ap+—T^ = ap + b, (8) 



P = t^P, 



T 



where w = (7 - 1), - 0.5 - 1.5 GeV , a = -(97^ - 247 + 16)/[9wo(7 
(3 = 0.718 and 7 ~ 1.318. For Eckert fluid, we insert Eq.[n]in Eq.Oto get 

2 In(iJ) - In {b-jH + 3aH^' 



t = 



-1 
3& 



27 arctan 



V 



^12q;6 



1)], 



(9) 



The results are depicted in left panel of Fig. [T] and compared with the non- viscous 
fluid, t = 1/2H. At b=0, i.e., ^ (x p, Eq.[9]turns to be 

- 3aH \nH + 3aH \n{3aH - 7) " 



t{H) = 2 



(10) 
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For IS fluid, let us first assume that the last term in rhs of Eq. [7] is too small 
compared to first term. With this assumption Eq. [5] can be solved and we get 
T cx R^/p, which has been discussed in Ref.^ as one of unphysical cases, in which 
T increases with the Universe expansion. Substituting Eq. [7]into Eq. [5] leads to 

H + |[1 + (1 - rh]HH + iH-{l + r)^ + %- 2)H^ + ^1h' = (n) 
z a HA 2 a 

which is a non-homogeneous and non-linear second type Abel differential equa- 
tion. If we would assume that the last two term in Ihs are vanishing, we then 
left with Bernoulli equation with the solutions H = and H w —2^/[3a{2 — 7)]. 
In both cases, H does't depend on t, i.e., static Universe. Again, if it is possible 
to express ^ depending on t, we get 2H + 3-/H'^ - 3C(t)H = 0, which is BernouUi 
equation and has the solution, H(t) — exp [| J S^{t)dt\ / {c + exp [^7 J £_{t)dt\ ) . The 
numerical solution for Eq. [Tl] is given in right panel of Fig. [TJ It strongly depends 
on the initial conditions for H and H. 



One of the direct implications of our assumptions is the relation t vs. T. From 
Eckert solution, p{t) ~ 3H^{t), and Eq. [8]we get 



t = 



\nT^^ - \n{-B -D + 2AT^)"-" - ln(-B + D + 2AT 

iCD 



■^2\B+D 



(12) 



where A = 3aP, B = 'S-f{a/3)-^/^ , C = 9T^, D = {B^ + 12AC)^^^. Eq. [TH is 
depicted in Fig. [2] and compared with the non-viscous fluid, t = cT^^. 




Fig. 1. Left panel: t vs. H{t) in Eckert fluid, Eq. |9] compared with the case when b in Eq. |8] 
vanishes and with non- viscous fluid, t = 1/2H. Right panel: numerical solution of H{t) in IS-fluid, 
Eq. 111! compared with non-viscous fluid. 
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Eckert Fluid 




Fig. 2. t vs. T in Eckert fluid, Eq. 1121 compared with the non-viscous fluid, t ocT 
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